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Radiative shocks (also called supercritical shocks) are high Mach number shock waves that pho- 
toionize the medium ahead of the shock front and give rise to a radiative precursor. They are 
generated in the laboratory using high-energy or high-power lasers and are frequently present in 
a wide range of astronomical objects. Their modelisation in one dimension has been the subject 
of numerous studies, but generalization to three dimensions is not straightforward. We calculate 
analyticaly the absorption of radiation in a grey uniform cylinder and show how it decreases with 
\R, the product of the opacity x an d of the cylinder radius R. Simple formulas, whose validity 
range increases when xR diminishes, are derived for the radiation field on the axis of symmetry. 
Numerical calculations in three dimensions of the radiative energy density, flux and pressure created 
by a stationary shock wave show how the radiation decreases whith R. Finally, the bidimensional 
structures of both the precursor and the radiation field are calculated with time-dependent ra- 
diation hydrodynamics numerical simulations and the influence of two-dimensional effects on the 
electron density, the temperature, the shock velocity and the shock geometry are exhibited. These 
simulations show how the radiative precursor shortens, cools and slows down when R is decreased. 
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I. INTRODUCTION 

In many astrophysical systems, the effects of radia- 
tion on hydrodynamics arc strong. This is the case 
with fast outflows and shocks, such as those encoun- 
tered in jets, bow shocks produced by the interaction of 
jets with the surrounding interstellar medium^ or radia- 
tive shocks arising in the envelope of pulsating evolved 
starsi 2 - Nowadays, these processes can be observed with 
an increasing level of details. For example, high angu- 
lar resolution imaging of jets produced by Young Stel- 
lar Objects^ shows complex structures of bright knots 
and shocks inside the jets. The improvement of obser- 
vation techniques, and more especially the combination 
of spectroscopy and high angular resolution techniques, 
will allow the study of the physical properties of shocks, 
in addition to the study of their complex shapes. This 
combination will be achieved by AMBERS on the Very 
Large Telescope Interferometer for instance. 



However, the interpretation of these new data should 
go together with the improvement and the development 
of new models or theories and new sophisticated nu- 
merical codes. Since these codes are complex, extensive 
testing made through code inter-comparisons and com- 
parisons with laboratory experiments are very useful.- 
Therefore, code benchmarking can be viewed as one mo- 
tivation for pursuing laboratory high-energy density ex- 
periments. These experiments can also be considered as 
the most relevant approach to check assumptions and to 
provide hints and new ideas to address open questions or 
issues such as radiative shocks in astrophysics £ 

High-energy density laboratory astrophysics 
(HEDLA) experiments are mostly driven on large- 
scale lasers&^&^ i 10 ' 11 or on Z-pinchesi 12 ' 13 In a first 
kind of experiment, one measures microscopic quantities 
required for the determination of the equation of state 
and the opacities of hot and dense matter found in 
giant planet interiors or in stellar atmospheres. In 
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a second kind of experiment, one seeks to reproduce 
astrophysical phenomena in the laboratory. The size of 
the experimental targets is of the order of a millime- 
ter, implying typical durations of a few nanoseconds, 
whereas the scale of astrophysical phenomena are 
between 15 and 25 orders of magnitude larger. However, 
scaling laws help to bridge laboratory experiments and 
astronomical phenomenal 14 ' 15 In order to be relevant 
for benchmarking, the appropriate energy requirements 
need to be satisfied and, more importantly, a set of 
accurate time and space-resolved diagnostics must be 
provided. Many experiments have been devoted to the 
study of typical radiative hydrodynamical situations, 
such as (i) radiative blast waves in the context of 
supernovae remnants , 16 ' 17 ' 18 (it) radiative precursor 
shocks waves ; 9 ' 10 ' 19 ' 20 with applications to the studies 
of stellar jets, pulsating stars , 21 ' 22 and accretion shock 
during star formation^ and (Hi) radiatively collapsing 
jets relevant for protostellar outflowsi^^ Reviews can 
be found in Refs. HEIll 

In this paper, we shall concentrate on the modeling 
of radiative-precursor shock wave experiments performed 
with high-power lasers. By radiative-precursor shock 
waves, we mean highly hypersonic shocks that are driven 
continuously by a piston. Because of the high shock ve- 
locity, the shocked medium is ionized and emits radia- 
tion, which in turn ionizes and heats the cold unshocked 
gas and leads to the apparition of a radiative precur- 
sor. The structure of the shock depends on the shock 
velocity^ 7 - For a given shock velocity, the temperature 
of the shocked matter and its mass density are solutions 
of generalized Rankine-Hugoniot equations— 8 - For weakly 
hypersonic shocks, the mass density varies continuously 
throughout the shock. At higher shock velocities, it be- 
comes discontinuous and then again continuous when the 
velocity increases because of the contribution of the ra- 
diation pressure and energy density. 27 ' 29 In the so-called 
supercritical regime, the mass density varies discontinu- 
ously, but the temperature remains constant through the 
shocks except in the shock front where the temperature 
spikes. The extent of this region is of the order of the 
photon mean free path. 

The strong coupling between radiation and hydrody- 
namics in these shocks is difficult to model because of the 
strong gradients and different length scales involved for 
hydrodynamics and radiation transport. Additionally, 
an accurate set of opacities and equations of state for a 
wide range of plasma conditions is needed. Departures 
from local thermodynamical equilibrium (LTE) are also 
expected to be important, especially in the shock front, 
requiring detailed calculations of the monochromatic ra- 
diation intensities. A complete and detailed study of the 
shock structure can therefore be achieved currently only 
in a restricted, ID geometry. However, the finite radial 
size R (in the direction perpendicular to the shock prop- 
agation) of actual shocks in the cosmos and in laboratory 
experiments may introduce departures from the ideal ID 
behavior. The purpose of this paper is to examine the 




FIG. 1: (Color online) Propagation of the shock (radius R) 
in the cell (radius R t ) filled with Xe. The laser beams are 
shown for illustrative purpose: during the experiment, the 
laser pulse is ended when the shock propagates in the cell. 



importance of these 3D geometrical effects. 

In a ID description, the radial extension R of the re- 
gion in the shock that emits radiation is considered to 
be infinite. The amount of radiation absorbed at one 
point of the configuration (in the precursor for example) 
is therefore overestimated by comparison with a descrip- 
tion where R is finite. This geometrical effect results in 
a "radiative energy loss" because a fraction of the energy 
is radiated radially and does not heat the radiative pre- 
cursor ahead of the shock. This effect will be dominant 
when the photon mean free path becomes large compared 
to R and may have a strong impact on the development 
of the radiative precursor and on the shape of the shock. 

Deficiencies in the current modeling of radiative shock 
experiments provide the motivation for a detailed investi- 
gation of the importance of these lateral radiative losses. 
For instance, Bouquet et alJ^ have recently reported on 
supercritical radiative shocks created with a high-power 
laser and the associated modeling work. While hydrody- 
namical simulations reproduce the main features of the 
experiment, questions remain regarding the precise un- 
derstanding of the formation and of the structure of the 
radiative precursor. 

In this paper, we report a study of the importance of 
radiative losses on the structure of shocks produced in 
laboratory experiments, aiming at understanding the de- 
ficiencies of ID models. Section |H] briefly recapitulates 
the results of Bouquet et alr^- and some of these deficien- 
cies. An analytical estimation of the effects of the finite 
size on the radiation distribution is given in Sec. IIII Al 
The numerical calculation of the spatial structure of the 
radiation field in a stationary case without coupling to 
the fluid is described in Sec. IIII Bl A full description of 
the radiative shock structure and of the radiation field is 
studied using the 2D Lagrangian radiative hydrodynam- 
ics code FCI,— with various boundary conditions for the 
radiation fSec. HV[) . We conclude (Sec.|Vj that it is essen- 
tial to account for multidimensional effects in modeling 
radiative shocks with finite radial extension. 
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II. RADIATIVE SHOCK EXPERIMENTS 

Several shock experiments have been performed with 
high energy density lasers to generate supercritical 
shocksi 9 ' 10 i 19 ' 20 i 32 In such experiments the laser energy 
is converted into kinetic energy of a piston, which drives 
a shock in a millimctric target of radius Rt (see Fig. [1]) 
filled with gas or low-density material. 

The extension R (shock radius) of the experimental 
shock in the direction perpendicular to its propagation 
along the x axis (see Fig. [1} is given by the size of the 
laser spot focused on the piston. The diameter 2R of the 
focal spot is equal to 500 /im [full width at half-maximum 
(FWHM)] with a 250 /im flat center where the intensity is 
constant. The shock has the same cylindrical symmetry 
than the laser focal spot. 

The target is impacted by the laser beams on one of its 
two ends consisting in a piston composed of several layers. 
The first layer is ablated during the 0.8 ns duration of 
the laser pulse. This ablation produces a shock wave 
that propagates in the x direction. The second layer is 
a radiative shield in titanium that prevents the radiation 
of the laser or the radiation emitted during the heating 
of the piston by the laser from penetrating the xenon gas. 
The shock finally breaks out in the tube filled with xenon 
gas at 0.1 or 0.2 bar. The shock velocities achieved in the 
gas are in the range 50 - 100 km s _1 and the duration of 
the experiment is 10 ns. 

The spatially averaged temperature of the shock is 
measured^ 3 , to be 10 - 20 cV and the electron density 3 ^ 
in the precursor is between 10 18 and 10 20 cm -3 . The 
precursor propagates ahead of the shock front with a ve- 
locity that equals up to twice the velocity of the shock 
front. 

The preparation and the initial analysis of the experi- 
ments were performed using one dimensional (ID) radi- 
ation hydrodynamics codes (FCI^ MULTI 3 ^) . While 
the experimental shocks deviate from ID geometry, these 
codes give an overall good agreement with the experi- 
ments (shock velocity, temperature of the shocked ma- 
terial). However, they are only moderately successful 
regarding the radiative precusor. They fail to reproduce 
the shape of the electron density profile in the precursor 
on the whole interval of values accessible to measure- 
ment (from 10 18 to 10 20 cm" 3 in Ref.Ql). The different 
codes provide a range of predicted velocities for the pre- 
cursor. For example, for a particular shot, the velocity 
at n e = 3 10 19 cm -3 is 300 km s _1 in the ID version of 
FCI and 120 km s" 1 in MULTI, the last one being in 
accordance with the measured velocity. It was reported 
by Vinci et alM- that the time evolution of the radial 
extension of the shock was well reproduced by the 2D 
version of FCI, therefore showing that 2D effects have to 
be considered. 

In the work presented here, we estimate the 2D effects 
and the consequences of the finite lateral size R of the 
shock on the radiation field. When investigating the ge- 
ometry of the problem, one needs actually to consider 
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FIG. 2: (Color online) Geometry of the configuration for the 
analytical model. The disk is located at x = 0. The limit 
angle (9 m i n = cos -1 /i m i n is the angle under which the disk is 
seen from the point M(x). 

the optical depth, t, which is the geometrical distance 
normalized by the photon mean free path. The opac- 
ity of the medium is therefore an essential parameter of 
the problem. Considering the overall agreement of the 
ID codes with the experimental results, we may assume 
that the current opacities are a reasonable approxima- 
tion. As a consequence, we will not examine that specific 
point further. 

III. RADIATION FIELD EMITTED BY A 
SHOCK 

A. Analytical estimation of the multidimensional 
effects 

We start with an analytical study of the radiation field 
generated by a disk in a grey and homogeneous cylinder 
with radius R and length L (Fig. [2]). Wc furthermore con- 
sider that the opacity x (cm -1 , the inverse of the photon 
mean free path) is uniform (independent of space) and 
we neglect scattering for sake of simplicity. The opti- 
cal depth in the x direction t(x) = J Q xdx, is the po- 
sition along the x axis normalized by the opacity. The 
cylinder length L is such that \L =100. The disk radi- 
ates an isotropic specific intensity I (J m s -1 sr" 1 ) 
and is characterized by its radius R or by its lateral 
optical depth tr = \R- The local source function S 
(J m~ 2 s _1 sr _1 ) is uniform in the cylinder and much 
lower than the specific intensity /„. It is zero outside of 
the cylinder. 

In this approach, wc would like to calculate the fraction 
of radiation energy emitted by the disk and received at 
any point M(x) of the cylinder axis x. In addition, it is 
aimed to examine the way this amount of energy varies 
by changing the disk radius. A brief presentation of this 
calculation has been presented in Ref . 1371 and an extensive 
version can be found in Ref. [2(]. We will use the three 
Eddington variables: the mean intensity J, the flux H, 
and the tensor K analogous to the pressure tensor. They 
are, respectively, related to the radiation energy density 
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E, the radiative flux F, and the radiative pressure tensor 
P through the relations (see Ref. [38h : 



./ 
H 
K 



(r) 




(1) 

where c is the speed of light in vacuum. In these expres- 
sions, the integration is performed over the solid angle u>, 
and J(r, n) is the specific intensity at position M(r) in 
the direction given by the unit direction vector n. Equa- 
tions |T]) arc frequency-dependent, but since we assume 
that the opacity \ is grey, wc will consider that the radia- 
tive quantities appearing in Eqs. ([1]) are integrated over 
frequencies. We will refer to J, H and K as the radia- 
tive moments since they correspond respectively to the 
zcroth-, first- and second-order moments of the intensity 
over angles. 

The formal solution^ of the radiation transfer equa- 
tion leads us to decompose the radiation field at position 
M(r) in two components: the contribution J D , H and 
K coming from the disk, and the emission Jg, Hg and 
Kg by the uniform medium. The total radiation field is 



(2) 



In this simple model, it is straightforward to calculate 
the radiation field existing at point M(r) on the x axis 
of the cylinder (r = 0). The contributions of the source 
function, Jg, H x $ and K xx s, are calculated by summing 
the emission coming from all the points of the cylinder. 
The contribution of the disk is given by 
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p e~ T ^dp, (3) 
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where we have considered only the x component of the 
flux and the xx component of the tensor K G . The quan- 
tity p is defined by p = cos 6, where 6 is the angle be- 
tween the direction n and the x axis. The integrals given 
by Eqs. ([3]) are exponential integrals in the case p mul = 
corresponding to a disk with a radius R going to infinity. 
For a finite radius of the disk, the specific intensity is in- 
tegrated over the solid angle Stt under which we see the 
disk from point M(r) (see Fig. [2]). 

The total value J(t) = J {t) + Jg(r) (and similarly 
for H and K) is plotted in Fig. [3] Because the general 
behavior of J, H x and K xx is similar, we can restrict 
our description mostly to the mean intensity J. The 
decrease of the radiative moments with tr enables us to 
understand the effect of the finite radial extension R of 
the medium on the structure of the radiation field. 

Far from the disk, the radiation emitted by the disk 
has been considerably absorbed, and the moments de- 
pend only on the local source function S. In the calcula- 
tions shown in Fig. [31 the ratio S/I Q = 10~ 4 . This value 
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FIG. 3: (Color online) Radiative moments (a) J/I , (b) H x /I 
and (c) K xx /I produced by a disk in a grey and homogeneous 
medium with S = 10 _4 7 o . The exact solution is shown for 
the infinite case (plain line with crosses) and for tr = 1, 1/10 
and 1/100 (stars, triangles and big dots). Also shown are 
the approximate solutions given by Eqs. <(Sj (dashed lines) 
and Eqs. ((5J times pmin (plain curves), which are given by 
Eqs. ((6} for J and H x . Positions of r eq (rij) are shown on 
frame (a) for tr = 1 and 1/100. 



corresponds to a system in which the radiation is due 
to Planck emission and the temperature of the medium 
is 10 times lower than the disk temperature. The flux 
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H x (t) goes to zero when r is greater than 1. An order 
of magnitude of Jg can be estimated by assuming the 
medium is a sphere of radius R, 

Js(r>l)~^/ dfi f * dr S e~ T ~ S (1 - e- TR ) . 
* J -i Ja 

(4) 

Though a rough estimate, it is quite correct in the cases 
considered here. For small Tr, the expansion of e~ TR 
gives Js(t ^> 1) ~ Str, an evaluation quite close to the 
exact value [see in Fig. [3] the value of J(r > 10) ~ Js ~ 
Str]. _ 

An interesting feature of the decrease of J with tr can 
be illustrated with Fig. [3j for tr smaller than one, it 
takes a distance of a few tr for J to be absorbed by a 
factor of 10, whereas, in the infinite case, where tr is 
greater than one, the radiation is absorbed tenfold after 
a propagation of a few mean free paths of photons, i.e. 
at t k 1 - 2. The absorption of the radiation field is 
therefore determined by the smallest of the characteristic 
scales. For the system studied here, the smallest scale is 
either the mean free path of photons l/x or R. 

In the case tr = 1 for which the two scales are the 
same, J is not far from the mean intensity Jooo calcu- 
lated in the infinite case. At position r = 2, J a is only a 
factor of four lower than Jo oo ■ Therefore, a characteristic 
lateral optical depth tr of a few is sufficient to consider 
the medium as optically thick and consequently nearly 
identical to a medium with an infinite lateral dimension. 

The other two moments for the disk radiation Hq x 
and K ox .x have similar profiles, and the three moments 
Jo, Hqx and K oxx have nearly the same value sufficiently 
far away from the disk, where it is possible to have an 
approximate solution of Eqs. At a distance from the 
disk of a few times tr, the solid angle <5f2 is small and 
Mmin « 1 - 2 /2 k 1 - (t r /t) 2 /2. The limiting value 
/i m i n is close to unity as soon as r is larger than 3 tr or 
4 tr and asymptotically goes to unity. 

The limits of the integrals in Eqs. ([3]) are 

(/■■) <t "- 0, ^.t(t) , '^(^) 

\ o xx ) \ A*min / 

(5) 

This asymptotic expression (represented in Fig. [3]) is 
useful to approximate the exact Eqs. © over a range 
[ T o j r oq] and to describe more simply the absorption of 
the radiative moments. Let us first estimate the lower 
value To- At position r = 5tr, the error on J D defined 
as the relative difference (J^pp™* - J° xact )/ J^ xact is, re- 
spectively, 6%, 0.6% and 0.06% for, respectively, tr = 1, 
0.1 and 0.01. We see on Fig. [3] that To is of the order 
of a few tr. Therefore, the approximation stands only 
far enough from the disk, where the mean intensity has 
been greatly reduced by the absorption. For example, 
J D (r = 5r fl )/J o (0) k 1/100 and 1/50 for respectively 
tr = 0.1 and 0.01. 

Let us define t c<1 (tr, I q / S) as the optical depth for 
which the disk radiation is equal to the local radiation. 



This is a straightforward way to separate the medium in a 
region r < r oq where the radiation field coming from the 
disk dominates and the region r > r cq where the disk 
does not have a direct influence on the radiation field. 
The transition between the two regions is sharp because 
the radiation field emitted by the disk decreases expo- 
nentially around r Gq . This r Gq can be easily calculated 
numerically. Its approximate position is shown Fig. [3Ja 
for tr = 1 and 1/100. We find that r eq decreases slower 
than tr since the local radiation Js decreases like tr 
when tr is small [see Eq. 0], whereas the disk radi- 
ation decreases like {tr/t) 2 exp(— r). We see in Fig. [3] 
that T cq ~ 5 for tr = 1 and r cq ~ 2 for tr = 1/100. 
Therefore, T cq decreases by less than a factor 5 when tr 
decreases by a factor 100. Then, in this particular cylin- 
drical configuration, the range [To,T eq ] over which the 
approximation ((5|) is good increases as tr decreases since 
To decreases proportionally to tr and r eq does not vary 
a lot. 

Equations (O show that the radiative moments on the 
x axis are absorbed by a factor e _T / /<mln and vary like 
(tr/t) 2 which, for a uniform opacity, is the solid angle 
(R/x) 2 under which the disk is seen from position r. 

Furthermore, we can calculate from Eqs. ([5]) the Ed- 
dington factor / = K oxx /J « Mmin an( i * ne anisotropy 
factor H ox /J « /-t m in- We conclude that the radiation 
field has a high degree of anisotropy at distances where 
the approximation is acceptable. 

The approximation ((5|) can be improved for J Q and 
H ox (but not for K oxx ) if they are multiplied by /i m i n , 




(6) 

The error for J a is then decreased by a factor of 5 to 
10, and this approximation is valid for a larger range 
than approximation ([5|) (see Fig. [3]), namely for values 
of t greater than approximately tr/2. For H OXl the 
improvement is very good since the error is then lower 
than tr for all values of r. It is therefore very useful for 
tr < 0.1, since the accuracy is then better than 10%. 
With this new approximation, the Eddington factor be- 
comes / = K xx/ Jo ~ Mmin and, when compared with 
the exact value, proves to be more accurate than /i^ lin . 



B. Numerical calculation of the radiation field 
generated by a planar shock 

1. Multidimensional radiative transfer with the short 
characteristics method 

The above analytical calculation is limited to the val- 
ues of the radiation field on the x axis in a uniform 
medium with a grey opacity. We now perform a multi- 
dimensional numerical computation of the radiation field 
generated at a given time by a planar shock. We vary 
the lateral size of the shock and study how the radiation 
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FIG. 4: (Color online) 3D and ID numerical computation 
of the flux H x (from the model presented in Sec. HUB 3|l for 
several values of R. The lines with symbols are the mean 
values, the dashed lines are the values at position r = 0, and 
the error bars are the mean deviations. The temperature T(x) 
of the shock is plotted. 



field is modified while all the other quantities are kept in- 
variant. The coupling of the radiation and the fluid will 
be studied in Sec. IIVI For the moment, we are only in- 
terested in the shape and the value of the radiation field 
in a nongrey description. 

The structure of the shock (temperature and density) 
along x is given by a one-dimensional radiation hydrody- 
namics code developed by J. P. Chieze and that has also 
been used in Ref. [l(J- It is chosen to be as close as possi- 
ble to the laboratory shocks we have observed^ in xenon 
(velocity ~ 60 km s _1 , P = 0.2 bar). The temperature 
is shown in Fig. |4j This ID shock structure is put in a 
three-dimensional (3D) grid by assuming a plane parallel 
shock structure. 

In the 3D radiative transfer calculation, the medium 
is assumed to be in LTE. The source function is there- 
fore assigned to be the Planck function. The opacity is 
approximated by a screened hydrogenic model, and cal- 
culated at n„ = 200 frequencies. 

The radiative transfer code calculates the specific in- 
tensity I(r, n) at each grid point of a 3D Cartesian grid 
(number of grid points: n x = 200, n y = 40, n z = 40) with 
the short characteristics method^ We use the cartesian 
coordinates because a solution of the radiation transport 
in these coordinates is easier to express and faster to 
compute than in the equivalent "natural" cylindrical sys- 
tem (see also Ref. |40| ). The smallest radial size of the 
shock for which the calculation has been performed is 
AY = AZ = 2R = 300 /mi, i.e., the size of the ex- 
perimental shock. For this value of R, and for a cell 
length Ax = 2 mm, the number of points for the angles 
is ?i M = 200, n v = 64 where the range of /j, is [—1, 1] and 
the range of ip is [0, 27r]. Solving the radiation transport 
in 3D is much more demanding than in ID not only be- 
cause of the larger number of grid points (n x , n y and 
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FIG. 5: (Color online) 2D plot of the mean in- 
tensity J(x,r) (J cm~ 2 s _1 ) for a planar shock of 
width R = 250 /jm. Isocontours are drawn for 
J = 2.10 s , 4.10 s , 6.10 8 , 8.10 s , 10 9 , 2. 10 9 ,4.10 9 and 6.10 9 
J cm - s . 



n z ) but also because of the required angular sampling 
(n v and n M ) which must be more refined as R decreases. 
Generally, this implies that some trade-off such as lim- 
iting the number of frequencies describing the opacity 
is necessary. In the ID case, 10 5 to 10 6 frequencies can 
be easily managed, whereas this number is drastically re- 
duced from 10 to 10 3 in 3D calculations (here, n„ = 200). 
In this calculation, the total number of mesh points is 
n x n y n z n v n v ~ 10 12 . 



2. Curvature of the radiation field 

The curvature of the radiation field is evidenced in 
Fig. [5] This curvature is due only to the finite lateral 
size of the medium since the source function is kept pla- 
nar, that is S(x, r) = B(T(x)), with B{T) the Planck 
function. The profiles of H x and K xx are similar to the 
profile of J. At position x = 500 pan in the radiative pre- 
cursor, at mid distance between the shock front and the 
ionization front, the value of the mean intensity on the 
boundary J(x,R) is roughly a factor of two lower than 
J{x,r = 0), the value on the axis. The medium is then 
in an intermediate state between optically thick and op- 
tically thin. At places x < 200 /im where the medium is 
optically thick, for example in the shock front, J is nearly 
constant, which makes it identical to the ID case. The 
medium becomes optically thinner as we go farther away 
from the shock front. The radiation field deviates from 
the planar geometry and less radiation is available to ion- 
ize the medium on the outer border of the ionization front 
(r ~ R) than on the axis. This leads the ionization front 
to also depart from a plane. This result is confirmed by 
the subsequent calculations of Sec. IIVI 
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3. Comparison between ID and 3D 

A comparison between the radiation field in 3D and in 
ID can be done by computing its value averaged on the 
surface AYAZ at each position x: 
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where 3 is either the mean intensity J, a component of 
the flux Hi or a component of the tensor , where i,j = 
x,y, z. 

The flux H x calculated with the 3D model for R = 150, 
250 and 500 (im and in ID for R = oo is shown Fig. 2] 
The error bars are standard deviations on the lateral sur- 
face and the dashed lines give the value at r = for 
each value of R. The magnitude of the radiation field de- 
creases with the lateral size of the shock, as was shown by 
our analytical approach. Although the radiation field is 
nearly constant in the ID infinite calculation, it decreases 
by one order of magnitude between the shock front (the 
spike of temperature) and the ionization front (where the 
temperature goes to zero) for R = 500 /im and by two 
orders of magnitude for R = 150 fxxa. We should note, 
however, that these calculations most likely overestimate 
the radiation field since the temperature profile is kept 
the same for all the values of R. But the decrease of the 
radiation field should result in a cooler and shorter pre- 
cursor, which will therefore emit even less radiation than 
the ID infinite- R profile. We are going to see in Sec. IIVI 
that when R decreases, the temperature in the precursor 
decreases too. 



IV. BIDIMENSIONAL TIME-DEPENDANT 
NUMERICAL SIMULATIONS 



In the general case, the determination of the radiation 
field [for example by solving Eqs. is difficult because 
the value of the optical depth r(r), which depends on the 
opacity x( r )> must be known. But the opacity depends 
on the local density, temperature and the ionization and 
excitation state of the medium, which can be modified 
by the radiation field. 

In the previous sections, we studied the consequences 
on the radiation field of the variation of the lateral size of 
the medium, without calculating the effect on the struc- 
ture of the medium. We now focus on the consequences 
of the bidimensional shape of the shock on the radia- 
tive precursor structure. We first compare ID and 2D 
numerical simulations with various boundary conditions 
with respect to radiation, and investigate the transition 
between these two kinds of calculations, stressing the role 
of the radial (transverse) radiative flux. In a second part, 
we study the 2D structure of the radiation field. 



Rt=500 transparent 
Rt=5()CI reflecting 
Rt=450 reflecting 
Kt=350 reflecting 
Rt=3(ltl reflecting 
Rt=250 reflecting 
ID 




■ Rt=500 transparent 
■Rt=500 reflecting 
■Rt=45Q reflecting 
■Rt=3J0 reflecting 

■ Rt=30C! reflecting 
■R[=250 reflecting 

■ ID 




-+■- 



FIG. 6: (Color online) Transition from the FCI-1D reference 
simulation to the FCI-2D reference case, using decreasing tar- 
get radius R t (in /im) and reflecting wall cells. Axis profiles 
of (a) electron temperature (eV), (b) electron density (cm -3 ) 
and (c) radiation energy density (J. cm" 3 ) are plotted at t = 
3 ns. 
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Mean Free Path (mm) 




FIG. 7: (Color online) 2D plot of Rosseland mean free path 
in mm, in the reference 2D case, at t = 5 ns. 



A. Transition from ID to 2D 

We perform time-dependant 2D radiation hydrody- 
namics simulations in cylindrical geometry with the code 
FCI<2I This code, devoted to the modeling of inertial 
confinement fusion, is therefore well suited to simulate 
the radiative shock experiments described in Sec. ITT1 from 
the laser pulse to the shock propagation in the Xe gas. 
The radial and temporal laser profiles can be well ap- 
proximated by constant functions, with Gaussian wings 
and 500 /mi FWHM. The value of the laser flux (3.10 13 
W/cm 2 at peak power in time) is chosen so that the 
velocity of the produced shock (55 km/s) matches the 
experiment shock velocity. 

We define a reference 2D simulation similar to a typical 
shot performed in the 2002 campaign. 10 i 32 The target 
diameter is 1 mm and the walls of the cell, which were 
made of quartz, are transparent to radiation. 

The transition between 2D and ID has been investi- 
gated by performing 2D simulations with walls reflecting 
radiation and decreasing the target radius Rt- However, 
the laser profile and the shock radius R are the same in 
all simulations. Reflective walls prevent lateral radiative 
losses at the boundary r = R t . A 2D simulation with 
reflective walls and R t equal to the laser spot radius R 
is nearly equivalent to a ID simulation. On the con- 
trary, a calculation with R t much larger than R clearly 
exhibits the 2D aspects of the problem. Increasing R t 
in these simulations is then comparable, although not 
strictly equivalent, to decreasing the shock radius in the 
model of Sec. lIIIl As stated earlier, the geometrical effect 
of decreasing R can be described as a radiative loss. 

This procedure emphasizes the importance of radia- 
tive losses both in simulations of radiative shocks and 
in related experiments by varying only one geometrical 
parameter. 




x (mm) 

FIG. 8: (Color online) Comparison of axis longitudinal radia- 
tive flux (J cm~ 2 s _1 ) with various boundary conditions, at t 
= 3 ns. 

In Fig. [6] a, we compare the axial temperature profiles 
at t=3 ns in the reference 2D simulation (lower curve) 
and the reference ID simulation (upper curve) computed 
with the ID version of the code^ Intermediate plots give 
electron temperature profiles for 2D simulations using re- 
flective walls and decreasing target radius Rt = 500, 450, 
350, 300, and 250 /im. Reflecting boundaries suppresse 
radiation losses and give higher radiative flux and energy 
in the precursor (Fig. [He), which then produce higher 
temperatures (Fig. [BJa) and electron densities (Fig. Ob). 
As the radius of the reflecting wall simulations decreases, 
the resulting temperature profile gets closer and closer 
to the reference ID profile. With reflecting walls and a 
350 fim target radius, the axial 2D electron density in 
the precursor is about half that of the ID value. With 
reflecting walls and a 250 /im radius, the 2D simulation 
yields an axial profile of the electron density quite simi- 
lar to the ID simulation, as expected. Radial flux losses 
thus seem to be the main source of differences between 
ID and 2D simulation results. 

One can also notice that the difference between the ref- 
erence 2D profiles and the others increases farther away 
from the shock front. This can be explained by the cumu- 
lative effect of radial losses on the radiative flux emitted 
by the shock front. 

A cartography of the Rosseland photon mean free path 
in xenon is shown Fig. [7] for the reference 2D simulation 
at t=5 ns, when the shock position is x ~ 250 /im. Only 
the xenon gas has been represented in the figures. The 
region in white color corresponds to the xenon gas, which 
has been swept up by the shock initially produced at 
x = 0. The radius of this region is about 300 /im at the 
origin x = whereas the radius of the laser focal spot is 
about 250 /im. We conclude, therefore, that the radius of 
the shock, R, increases with time while it is propagating 
to the right in the xenon still at rest. In addition, the 
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FIG. 9: (Color online) 2D profile of the electron temperature 
in the gas in the reference 2D case, at t = 5 ns. 

curvature of the shock front is clearly evidenced. 

The values of the Rosseland mean free path 1/x in 
the precursor (400 /im < x < 800 /mi) range typically 
between 200 /im and 1 mm (or less), which correspond, 
respectively, to tr = \R ~ 1-25 and 0.25, whereas they 
strongly decrease in the shocked xenon and the shock 
front (250 /im < x < 400 /im), where t r « 2.5 for \ « 0.1 
mm. With these characteristic values of tr, the shock 
might be considered as nearly ID in the shock front and 
3D in the precursor. 

The effects of boundary conditions on the radiation be- 
come less important as the radius of the target increases, 
as seen in Fig. [5J The longitudinal flux calculated for 
transparent walls increases as the radius of the target 
increases, and the volume of gas ionized by radiation in- 
creases too. The Rosseland photon mean free path (see 
Fig. [7]) lies typically between 200 /im and 1 mm or less in 
the precursor. Therefore, most of the radiation emitted 
by the shock has been absorbed by the gas before reach- 
ing the walls of the cell. The radiation emitted by all 
the ionized gas contributes to the longitudinal flux which 
increases until the maximum volume of ionized gas is ob- 
tained in the precursor, ahead of the shock. 

On the other hand, the longitudinal flux calculated 
with reflective boundaries decreases as i?t increases be- 
cause the radiation that reaches the walls when the radius 
is small enough is reflected and then contributes again to 
the flux and to the ionization of the gas. The limit sit- 
uation is reached when the boundary is sufficiently far 
away so that most photons are absorbed in the gas. 

Another consequence of the smaller amount of radia- 
tion emitted in 2D than in ID is the smaller precursor 
velocity. The precursor velocity, defined as the speed of 
a chosen isocontour value of electron density, decreases 
from 350 km s _1 in the ID simulation to 100 km s _1 
in the 2D reference case, for n e = 10 19 cm -3 . The ID 
simulations thus also overestimate precursor velocities in 
a systematic way. 
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FIG. 10: (Color online) 2D plot of the (a) longitudinal and 
(b) radial fluxes (in J cm' 2 s' 1 ), and (c) their ratio in the 
reference 2D case, at t = 5 ns. 



B. Bidimensional structure of the radiation field 

Figure [5] shows that the electron temperature is much 
higher in the shock front than in the precursor: 25 eV 
versus 5 eV and less. Assuming that emissivity scales 
as T 4 like a black body, this makes precursor emissivity 
mostly negligible compared to the emissivity of the shock 
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front. Since in addition the average mean free path in 
the precursor is similar to the cell and precursor length 
scales, we can expect the radiation field to be effectively 
dominated by the radiation coming from the shock front. 
This will not be true anymore far enough from the shock 
front. The analysis of Fig. [8] done in Sec. IIV Al shows 
that, for i? t > 500 /im, the radiation from the volume of 
the gas (the precursor) becomes important at distances 
from the shock front greater than about 1 mm. 

Let us consider the radial and longitudinal components 
of the radiative flux, plotted in Fig. [TUla and Fig. [TT31b. 
On the x axis, the longitudinal flux is maximum, whereas 
the radial flux is zero because of the cylindrical symme- 
try. The radial flux is most important near the shock 
front. On Fig.Qjllc showing the ratio between both fluxes, 
we see that the medium can be separated in two regions 
with different regimes. On the dashed line, the two com- 
ponents are roughly equal, while the radial flux domi- 
nates above this line and the longitudinal flux dominates 
closer to the axis, below this separation. 

The radiation flux field appears, therefore, very 
anisotropic near the axis. It is strongly oriented in the 
shock propagation direction. This structure of the radi- 
ation field is very similar to that obtained in the model 
presented in Sec. IIIII which had a planar shock struc- 
ture and transparent walls. This suggests that the 2D 
radiation at r < 250 /im could be well modelled by a 
planar uniform source with finite radius R (representing 
the shock front) and transparent walls. 

V. CONCLUSION 

Geometrical effects on the radiation field are impor- 
tant when the radial optical depth = \R is lower 
than unity. Analytical arguments explain the decrease 



of the energy density, flux, and pressure when the radius 
R of the emitting surface diminishes. They show that, 
when R is lower than the photon mean free path 1/x, 
the radiation is absorbed roughly by a factor of ten at 
a distance x from the emitting surface of a few R. The 
dimension of the heated region therefore scales like R in 
this simple model. Approximate solutions of the radia- 
tion energy density, flux, and pressure are found in the 
asymptotic limit x > R. Further studies are needed to 
estimate their relevance in more complexe models. Veri- 
fications of the radiative transfer code can also be made 
by comparison with the analytical solution. 

Numerical calculations of radiative shocks confirm that 
the radius of the shock must be taken into account when 
tr = X-R < 1 even though % is a mean opacity (the 
Rosscland mean in the present case), which is a very 
rough representation of the actual frequency dependence 
in the radiative transfer calculations. One-dimensional 
calculations for such a complex system overestimate the 
amount of energy radiated and transfered ahead of the 
shock. This results in overestimating the energy depo- 
sition in the precursor. Therefore, all the properties of 
the radiative precursor are overestimated: its velocity by 
a factor 3, the temperature and electron density by one 
order of magnitude, and the extension by a factor that 
depends on the electron density, but can be up to an 
order of magnitude. These values are for the properties 
of the precursor far enough from the shock front, say at 
more than a photon mean free path. Moreover, the shape 
of the radiation field in the precursor and the structure 
of the temperature and electron density in the precur- 
sor depart from a plane parallel geometry. Our work, 
therefore, emphasizes the need of considering 3D or 2D- 
cylindrical numerical simulations for modeling radiative 
shocks in the laboratory and in cosmic settings. 



* Electronic address: |sleygnac@cp . dia s .ie | 

1 B. Reipurth, S. Heathcote, J. Morse, P. Hartigan, and 
J. Bally, Astron. J. 123, 362 (2002). 

2 E. Huguet and J.-P. J. Lafon, Astron. & Astrophys. 324, 
1046 (1997). 

3 C. S. Gil, E. M. Thiebaut, P. Garcia, and M. Schoeller, 
in New Frontiers in Stellar Interferometry, Proceedings of 
SPIE, edited by Wesley A. Traub. Belhngham, WA: The 
International Society for Optical Engineering (2004), vol. 
5491, pp. 1742-1747. 

4 A. C. Calder, B. Fryxell, T. Plewa, R. Rosner, L. J. Dursi, 
V. G. Weirs, T. Dupont, H. F. Robey, J. O. Kane, B. A. 
Remington, et al., Astrophys. J. Suppl. Ser. 143, 201 
(2002). 

5 B. T. Draine and C. F. McKee, Ann. Rev. Astron. Astro- 
phys. 31, 373 (1993). 

6 J. Grun, J. Stamper, C. Manka, J. Resnick, R. Burris, 
J. Crawford, and B. H. Ripin, Physical Review Letters 66, 
2738 (1991). 

7 D. R. Farley, K. G. Estabrook, S. G. Glendinning, S. H. 



Glenzer, B. A. Remington, K. Shigemori, J. M. Stone, R. J. 
Wallace, G. B. Zimmerman, and J. A. Harte, Physical Re- 
view Letters 83, 1982 (1999). 

8 K. Shigemori, R. Kodama, D. R. Farley, T. Koase, K. G. 
Estabrook, B. A. Remington, D. D. Ryutov, Y. Ochi, 
H. Azechi, J. Stone, et al., Phys. Rev. E 62, 8838 (2000). 

9 J. C. Bozier, G. Thiell, J. P. Le Breton, S. Azra, M. De- 
croisette, and D. Schirmann, Phys. Rev. Lett. 57, 1304 
(1986). 

10 S. Bouquet, C. Stehle, M. Koenig, J.-P. Chieze, A. Benuzzi- 
Mounaix, D. Batani, S. Leygnac, X. Fleury, H. Merdji, 
C. Michaut, et al., Phys. Rev. Lett. 92, 225001 (2004). 

11 N. C. Woolsey, Y. Abou Ali, R. G. Evans, R. A. D. 
Grundy, S. J. Pestehe, P. G. Carolan, N. J. Conway, R. O. 
Dendy, P. Helander, K. G. McClements, et al., Physics of 
Plasmas 8, 2439 (2001). 

12 J. E. Bailey, G. A. Chandler, D. Cohen, M. E. Cuneo, 
M. E. Foord, R. F. Heeter, D. Jobe, P. W. Lake, J. J. 
Macfarlane, T. J. Nash, et al., Physics of Plasmas 9, 2186 
(2002). 



11 



13 S. V. Lebedev, J. P. Chittenden, F. N. Beg, S. N. Bland, 
A. Ciardi, D. Ampleford, S. Hughes, M. G. Haines, 

A. Frank, E. G. Blackman, et al., Astrophys. J. 564, 113 
(2002). 

14 D. D. Ryutov, R. P. Drake, J. Kane, E. Liang, B. A. Rem- 
ington, and W. M. Wood-Vasey, Astrophys. J. 518, 821 
(1999). 

15 D. D. Ryutov, B. A. Remington, H. F. Robey, and R. P. 
Drake, Phys. Plasmas 8, 1804 (2001). 

16 K. A. Keilty, E. P. Liang, T. Ditmire, B. A. Remington, 
K. Shigemori, and A. M. Rubenchik, Astrophys. J. 538, 
645 (2000). 

17 J. M. Laming and J. Grun, Phys. Rev. Lett. 89, 125002 
(2002). 

18 M. J. Edwards, A. J. MacKinnon, J. Zweiback, K. Shige- 
mori, D. Ryutov, A. M. Rubenchik, K. A. Keilty, E. Liang, 

B. A. Remington, and T. Ditmire, Phys. Rev. Lett. 87, 
085004 (2001). 

19 P. A. Keiter, R. P. Drake, T. S. Perry, H. F. Robey, B. A. 
Remington, C. A. Iglesias, R. J. Wallace, and J. Knauer, 
Phys. Rev. Lett. 89, 165003 (2002). 

20 S. Leygnac, Ph.D. thesis, Universite Paris XI (2004). 

21 Y. A. Fadeyev and D. Gillet, Astron. & Astrophys. 333, 
687 (1998). 

22 Y. A. Fadeyev and D. Gillet, Astron. & Astrophys. 354, 
349 (2000). 

23 C. Stehle and J.-P. Chieze, in SF2A-2002: Semaine 
de I'Astrophysique Francaise, Paris, France, edited by 
F. Combes and D. Barret (EDP Sciences, Les Ulis, 2002), 
p. 493. 

24 B. A. Remington, D. Arnett, R. P. Drake, and H. Takabe, 
Science 284, 1488 (1999). 

25 H. Takabe, Progress of Theoretical Physics Supplement 
143, 202 (2001). 

26 B. A. Remington, Plasma Physics and Controlled Fusion 
47, A191 (2005). 

27 S. Bouquet, R. Teyssier, and J.-P. Chieze, Astrophys. J. 
Suppl. Ser. 127, 245 (2000). 

28 C. Michaut, C. Stehle, S. Leygnac, T. Lanz, and 
L. Boireau, European Physical Journal D 28, 381 (2004). 

29 C. Michaut, C. Stehle, L. Boireau, and S. Leygnac, in Third 



international conference on Inertial Fusion Sciences and 
Applications, Monterey, Ca., USA, edited by B. Hammel, 
D. Meyerhofer, J. Meyer-ter-Vehn, and H. Azechi (Ameri- 
can Nuclear Society, La Grange Park, 2003), pp. 954-957. 

30 Y. B. Zel'dovich and Y. P. Raizer, Physics of shock waves 
and high-temperature hydrodynamic phenomena (Mineola, 
New- York : Dover Publications Inc., edited by Hayes, Wal- 
lace, D. and Probstein, Ronald, F., 2002). 

31 FCI1 and FCI2 are respectively ID and 2D codes devel- 
opped at the Commissariat a l'Energie Atomique (CEA, 
France). See more details in Dautray R. & Watteau J.- 
P., La Fusion Thermonucleaire par Laser (Eyrolles, Paris, 
1993). 

32 X. Fleury, S. Bouquet, C. Stehle, M. Koenig, D. Batani, 
A. Benuzzi-Mounaix, J.-P. Chieze, N. Grandjouan, J. Gre- 
nier, T. Hall, et al., Laser and Particle Beams 20, 263 
(2002). 

33 T. Vinci, M. Koenig, A. Benuzzi-Mounaix, C. Michaut, 
L. Boireau, S. Leygnac, S. Bouquet, O. Peyrusse, and 
D. Batani, Physics of Plasmas 13, 010702 (2006). 

34 S. Leygnac, C. Michaut, C. Stehle, T. Vinci, M. Koenig, 
A. Benuzzi-Mounaix, S. Bouquet, F. Thais, and H. Merdji, 
in Third international conference on Inertial Fusion Sci- 
ences and Applications, Monterey, Ca., USA, edited 
by B. Hammel, D. Meyerhofer, J. Meyer-ter-Vehn, and 
H. Azechi (American Nuclear Society, La Grange Park, 
2003), pp. 958-961. 

35 R. Ramis, R. Schmalz, and J. Meyer- Ter-Vehn (1986), rep. 
MPQ110, Max-Planck Institut fur Quantenoptik, Garch- 
ing, Germany. 

36 R. Ramis, R. Schmalz, and J. Meyer- Ter-Vehn, Comp. 
Phys. Comm. 49, 475 (1988). 

37 S. Leygnac, L. Boireau, C. Michaut, T. Lanz, C. Stehle, 
and S. Bouquet, Journal de Physique IV 133, 453 (2006). 

38 D. Mihalas, Stellar atmospheres (San Francisco, 
W. H. Freeman and Co., 1978), 2nd ed. 

39 P. Kunasz and L. H. Auer, Journal of Quantitative Spec- 
troscopy and Radiative Transfer 39, 67 (1988). 

40 M. van Noort, I. Hubeny, and T. Lanz, Astrophys. J. 568, 
1066 (2002). 



